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S3 ■ I. INTRODUCTION 

^ There is a famous theorem named in honor of Arthur Cayley and William Hamilton in linear algebra, which asserts 

, that any n x n matrix A is a solution of its associated characteristic polynomial xa Q • In the popular SU (2) x SU{2) 
and SU{3) x SU{3) chiral perturbation theories, the Cayley-Hamilton theorem has been used to eliminate redundant 
terms and keep only independent pieces of the chiral Lagrangian 0, 0, 0j Hi • The trace formulae of the Cayley-Hamilton 
J> ■ theorem for 2x2 and 3x3 matrices are respectively 

- tr(A)A + i[tT2(A) - tr(A2)] = , (1) 

O ■ - tr(A)A2 -t- i[tr2(A) - tr(A2)]A - -[tT^{A) - 3tr(A2)tr(A) + 2tr(A3)] = . (2) 

■ 2 6 

And a recursive algorithm to compute the characteristic polynomial coefficients of a generic nx n matrix as functions 
. of the traces of its successive powers was given by Silva nearly ten years ago [1] . The work of Ref . obtained the 
recurrence relations for the coefficient sequence by using an improved version of the remainder theorem and some 
Qh! additional results. 

^ • In this paper, with the knowledge of combinatorics, we give a new proof for these recurrence relations, and then we 
further obtain an explicit formula for the generic term of the coefficient sequence, which yields the trace formulae of 
Cayley-Hamilton's theorem with all coefficients explicitly given, and which implies a byproduct, a complete expression 
for the determinant of any finite-dimensional matrix in terms of traces of its successive powers. We also discuss some 
applications to chiral perturbation theory and general relativity. 
' This paper is organized as follows. In Sec. [Tll by using combinatorics, a new proof for the recurrence relations of 
the characteristic polynomial coefficients is given and the explicit expression for the generic term of the coefficient 
sequence is obtained, which improves the results of previous works. In Sec. IIIIl a detailed discussion of the applications 
of the Cayley-Hamilton theorem to n-flavor chiral perturbation theories up to n = 5 is presented. In Sec. llVi as a 
further application, the determinant of the metric tensor in 10 dimensional spacetime is explicitly expressed by the 
traces of its successive powers. Sec. |V]is devoted to conclusions. A check for the Cayley-Hamilton relations using an 
alternative approach is presented in App. \^ while a PASCAL code for the computation of the generic term of the 
coefficient sequence is given in App. [BJ 



X 



II. TRACE FORMULAE 



Let A be any complex nxn matrix with eigenvalues Ai, A2, ■ • ■ , A„ (possibly complex and identical). Then its trace is 
tr(A) = Ai-|-A2-I-- • --I- A„, and the characteristic polynomial is Xa(A) = det(AI— A) = (A— Ai)(A— A2) • • • (A— A„), where 
I is the nxn unit matrix. And A being similar to a Jordan normal form implies that tr(A™) = A™ + A™ + • • • + A™ 
for any positive integer m. Following the notation of Ref. the traces are denoted by Tm = tr(A"') here and 
henceforth. Now we are ready to give a new proof for the theorem relating the coefficients of Xa(A) to the traces of 
A™ (m = 1,2, • • • 
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Theorem 1 //xa(A) = A" + -DiA" ^ + ■ • • + D„_iA + Z?„ is the characteristic polynomial of the complex n x n matrix 
A and T„, = tr(A'"), then 

mAn + An-iTi + An-2T2 + • • • + + T,„ = , m=l,2,---,n. (3) 

Proof. Let Ai, A2, • • • , A„ be the n eigenvalues of the matrix A. Then xa(A) = A" + £>! A"^^ + • • • + _D„_i A + Z).,i = 
(A-Ai)(A-A2)---(A- A„) yields 

Dr,i^ (J m , crm= ^ Ai^ Ai^ • • • Ai^ , m = 1, 2, ■ ■ ■ , n . 

il<i2<---<im 

Multiplying CTm^i by (Ai + A2 + • • • + A„), we have 

<^m-lTl = (Ai + A2 + • • • + A„) • -^41 ^^12 ' ' ' -^iin-l 

^l<^2<•••<^77^-l 
n n 

i=l jl<h<---<'jm-2 i=l jl<j2<---<jm-l 

Now let us show that the second term of the last line of the above equation contains m times of am ■ For each specific 
term of am, say A1A2 ■ ■ ■ Am, it is included m times in the second term of the last line of the above equation, i.e., it 
is contained by the m terms as follows 

Xi ■ Ajj Aj2 • • • , z = 1, 2, • • • , m . 

jl<j2<---<jm-l 



Thus, 



am-lTl - '^Xi ■ X! • • • 



i=l jl<j2<---<jm~2 



Next, multiplying am-2 by (Af + A2 + h A^) gives 

n n 



Likewise, 



i=l jl<j2<---<im-3 *=1 Jl<j2<---<j,Ti-2 



i=l il<j2<---<Jm-4 i=l jl<j2<---<jm-3 



a,Tm-, = E^r+E^r'E^^ = ^™+E^r'E^^- 

i—1 i—l j^i i—1 j^i 

Combing the above equations yields 

fm-iTi — am-2T2 + (^m-^T^ ~ am-iTi + • • • + {—)™aiTm-i = mam + {—)'^Tm , 

which results in 

mDm + Dm-iTi + Dm-2T2 + ■■■ + D^Tm-i + r,„ = . □ 
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From Theorem [U the characteristic polynomials for n x n matrices can be recursively obtained as follows 
n=l , Xa(A) = a - Ti , 
n = 2 , 



1 



Xa(A) = A3 - TiA^ + i(T2 - T2)A - ^(Tf - STiTa + 2r3) , 



(4) 



It is interesting to note that the characteristic polynomial coefficients Dm (I < m < n) as polynomial functions of the 
traces Tj^ (k = 1,2, . . . ,n) are formally unchanged when n increases, though are obviously different for distinct n. 
As a consequence, we can regard Dm (m = 1,2, . . .) as an infinite sequence. Theorem [1] gives 



2D2 = -{T1D1+T2) , 

3D3 = -{T1D2 + T2D1+T3) , 



(5) 



which leads to 



Di + 2xD2 + Sx^Da + • • • 

-Ti - x{TiDi + T2) - x^{TiD2 + T2D1 + T3) 

Ti{l + xDi + x^D2 + ---) + xT2{l + xDi+ x^D2 + ■■■) 



(6) 



If we define the generating functions /(.t) and g{x) for the infinite sequences {Dn} and {Tk+i} respectively as follows 



f{x) = J2 ^"^n (^itli ^0 = 1), g{x) = ^"Tk+i , 



Tl = 



k=0 



(7) 



then 



T 



1 



jn— 1 



x=0 



(n - 1)! dx"- 



(" = 1,2,...) 



x=0 



and Eq. ^ implies 



or 



which gives 



fix) = cxp 



l{x)fix) , 



g{t)At 



f{x) = exp 



^ n 

n=l 



n 



n=l 



00 00 



HE 

n— 1 p„— 



=0 (pi,P2,...,p„)e5„ ™=i 



T 
m 



(8) 



(9) 



(10) 



(11) 



where the set 5„ is defined by all nonncgative-integer solutions {{pi,P2, ■ ■ ■ , Pn)} of the equation pi+2p2+- • ■+npn = n. 
Comparing Eq. (jlip with Eq. ([7]), we obtain the expression for a generic term of the sequence {Dn} as follows 



1 

(pi,P2,---,Pn)65„ "1=1 



T 
m 



(12) 
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Example 1. Consider the computation of D^. All nonnegative integer solutions {pi,P2iP3jP4) for the equation 
Pi + 2p2 + 3p3 + 4p4 = 4 and their corresponding terms are shown in Table [H Adding up all these terms gives 

Ti T1T3 T?T2 T? Tf , , 



TABLE I: All nonnegative integer solutions (pi,p2,P3,P4) for the equation pi + 2p2 + 3p3 + 4p4 = 4. For each solution, there 
is a corresponding term -^{—T-iY''--^{—'^Y'^-^{ — ^Y'^-^{ — ^Y*. Summing of all the corresponding terms yields Dn. 



Pi 


P2 


Pz 


Vi 


corresponding terms 











1 




1 





1 
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1 
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Now we state without proof the well-known theorem in linear algebra: 

Theorem 2 (Cayley-Hamilton) Any n x n matrix A obeys its own characteristic equation, that is, xa(A) = 0. 

Combining Eq. (|12p with Theorems [l] and [2l we obtain the following improved results: 

Proposition 1 For any complex n x n matrix A, its characteristic polynomial is XaW = A" + DiX"^^ + • • • H 
Dn-i^ + Dn, and xa(A) = A" + DiA"^^ + • • • + D„-iA. + D„ = where the coefficients are given by 



E n 

(pi,P2,---,Pfc)e5fc m=l 



1 



P71 



T 

TO 



(fc = 1,2,--- ,n) 



(14) 



with standing for tr(A"') and Sk being a set including all nonnegative-integer solutions {{pi,P2, ■ ■ ■ iPk)} of the 
equation pi + 2p2 + • • • + kpk = k. 

Note that the determinant of a matrix A is the product of all the eigenvalues A„i (to = l,2,...,n), that is, 
det(A) = (T„ = 0™=! Am- Thus, Proposition [T] implies the following byproduct results: 

Proposition 2 For any complex n x n matrix A, its determinant is given by 



det(A) = (-)"D„ = (-)" 



1 

(pi,P2,---,Pn)e5„ m=l 



T 

TO 



(15) 



with Tm standing for tr(A™) and Sn being a set including all nonnegative-integer solutions {{pi,P2, ■ ■ ■ ,Pn)} of the 
equation pi + 2p2 + • • • + np„ — n. 

III. APPLICATIONS TO CHIRAL PERTURBATION THEORY 

The trace formulae of the Cayley-Hamilton theorem can be used to build relations between traces of finite dimen- 
sional matrices. Let us study, case by case, the trace relations ofnxn matrices for n = 2,3, 4, and 5. In the following, 
(A) is employed to denote the trace of a matrix A. 



1. The n = 2 case 
For any 2x2 matrix A, Cayley-Hamilton theorem reads 

= xa(A) = A^ - (A)A + i((A)2 - (A^)) 



(16) 



Multiplying Eq. ^TE\i by A and taking the trace results in 

= F(A)^(A3)-^(A2)(A) + i(A)^ (1 

Substituting A = AiAi + A2A2, where A1.2 are arbitrary parameters, into Eq. (jl7[) one finds 

= F(AiAi + A2A2) 

= A?i^(Ai) + Aii^(A2) + A?A2Fi2(Ai, A2) + A2AiF2i(Ai, A2) 

= A?A2i^i2(Ai,A2)+A2AiF2i(Ai,A2) . (1 
Now, substituting A = AiAi + A2A2 + A3A3 into Eq. (fTT]) leads to 
= i^(AiAi + A2A2 + A3A3) 



= J2^fP{^^) + J2 >H^iF,j{M,A,) + \]KFj,{A,,Aj) +AiA2A3Fi23(Ai,A2,A3) 

2—1 i<ij 

AiA2A3Fi23(Ai, A2,A3) , (1 

where the last equality comes from Eqs. (fTT)) and ITSl) . Thus. Fi23(Ai, A2, A3) = 0, that is. 

0= (A1A2A3)- (AiA2)(A3) + (Ai)(A2)(A3) . (2 

2 perm 3 perm 

Likewise, multiplying Eq. (fT6|) by and subsequently taking the trace gives 

= (A-*) - (A^^)(A) - \{AY + ^(A2)(A)2 . (2 
Then inserting A = Yll^i ^i-^i i^i^o Eq. (pij) and comparing the coefficients of A1A2A3A4, one finds 

0= (AiA2A3A4)-J Y (AiA2A3)(A4)- Y (A1A2) (A3A4) + ^ Y (AiA2)(A3)(A4) . (2 

6 perm 8 perm 3 perm 6 perm 

In particular, if (A.;) = 0, Eq. ([22|) becomes 

= ^ (A1A2A3A4) - Y (AiA2)(A3A4) , (for (A,) = 0) , (2 

6 perm 3 perm 

whose explicit form is 

(A1A2A3A4) + (A1A3A2A4) + (A2A1A3A4) + (A2A3A1A4) + (A3A1A2A4) + (A3A2A1A4) 
= (A1A2) (A3A4) + (Ai A3) (A2A4) + (Ai A4) (A2A3) , (for (A,) = 0) . (2 

Then further taking Ai = A2 = A and A3 = A4 = B in Eq. ([M]) . one obtains 

4(A2b2) + 2(ABAB) = (A^) (B^) + 2(ABf , (for (A) = (B) = 0) . (2 

Eventually, multiplying Eq. (|16p by A'^ and the same procedure as above results in 

0= Y (A1A2A3A4A5A6) - ^ Y (AiA2A3A4A5)(A6)-^ Y (AiA2A3A4)(A5A6) 

120 perm 144 perm 90 perm 

+ ^ Y (AiA2A3A4)(A5)(A6) . (2 

90 perm 

Eqs. ((20|) . (|22p and ((26)) are useful relations for the construction of any 2-flavor chiral perturbation theory. 

Example 2. Consider a 2-flavor chiral perturbation theory with the chiral (global) symmetry breaking patte 
[SU{2)l X SU{2)r X U{\)v]/[SU{2)v x U(\)v]- The basic building block of the theory is usually chosen to be 
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unitary unimodular 2x2 matrix field U, which transforms as f/ RUL^ under the chiral symmetry group with 
R e SU{2)ji and L e SU{2)l. If we define = W{D^U) which transforms Uke -> LX^.L'^. As det C/ = 1, then 

= a^detJ/ = a^(e'°sdctc/) ^ ei°s^'=*^a^(logdet [/) = ^^(trlogC/) = iY{U^d^U) = ti{U^D^U) = tr(X^) . (27) 

In the chiral Lagrangian, the chiral symmetric terms that consist of four X^'s arc {X^X^^X^X'-') , {Xf^^X^X^^X'^) , 
(X^X'^)'^, and (X^^X^)^. ExpHcitly, the first two terms are 



(X^X^'X^X"^) = J2 gUX'X'X'X') + 9^^9kk ( (X^X'X^X^) + (X'^X'^X'X') 

i i<k ^ 

= 5]gf,(X'X'X'X') + 2^5,,gfcfc(X*X'X'=X^) , (28) 

i i<k 

(X^X.X^X"^) ^ ^g2^X^^^X'X*)+2^5,,gfcfc(X*X'=X'X'^) . (29) 



i<k 



Note that there is no summation on repeated indices i and k unless specially stated by a summation symbol, 
Then, a linear combination of these two terms gives 

2(X^X^X^X'') + (X^X^X^X"^) = 3^.g|(X*X*X*X*) +Y9i^9kk ( 4(X*X*X'^'X'^') + 2(X*X'=X*X'=) j 

i i<k ^ ^ 

^ ^ g2 j^.)2 ^ ^ ("(X^X*) (X'=X'=) + 2(X*X'=) (X'X'^)') 

1 



-(X^X^)2 + (X^X,)2, (30) 

where the second equality comes from Eqs. ([?!]) and ([^5]) . From Eq. one finds that at least one of the four 4-Xp 

terms can be determined by the other three terms. A check for Eq. (j30p by an alternative method is given in App. 

m 

2. The n — 3 case 
For any 3x3 matrix A, Cayley-Hamilton theorem reads 

= xa(A) = A^ - (A)A2 + i((A)2 - (a2))A - 1((A)3 - 3(A2)(A) + 2(A3)) , (31) 
from which, the same procedure as the n = 2 case yields 

= ^ (A1A2A3A4)- Y (AiA2A3)(A4) - Y (AiA2)(A3A4)+ Y (A1A2) (A3) (A4) 

6 perm 8 perm 3 perm 6 perm 

-(Ai)(A2)(A3)(A4) , (32) 
0= Y (A1A2A3A4A5A6) - ^ Y (AiA2A3A4A5)(A6)-^ Y (AiA2A3A4)(A5A6) 

120 perm 144 perm 90 perm 

+ ^ Y (AiA2A3A4)(A5)(A6) - Y (AlA2A3)(A4A5A6) 

90 perm 40 perm 

+ i Y (AiA2A3)(A4A5)(A6)-i Y (AiA2A3)(A4)(A5)(A6) . (33) 

120 perm 40 perm 

Eqs. ([311) and ([551) ^-^'^ useful relations for the construction of any 3-flavor chiral perturbation theory, such as 
[SU{3)l X SU{3)r X U{l)v]/[SU{3)v x U{l)v] and 5'C/(3)/5'0(3) chiral perturbation theories. 



3. The n = A case 
For any 4x4 matrix A, Cayley-Hamilton theorem reads 

= xa(A) = A4-(A)A-vi((A)2-(A2))A2-i((A)-^-3(A2)(A)+2(A3))A 

+ (^-i(A^) + i(A-^)(A) + i(A2)2-i(A2)(A)2 + l(A)4^ , (34) 
from which, the same procedure as the n = 2 case yields 

0= (A1A2A3A4A5A6) - ^ (AiA2A3A4A5)(A6) - (AiA2A3A4)(A5A6) 

120 perm 144 perm 90 perm 

+ ^ E (AiA2A3A4)(A5)(A6)- Y (AiA2A3)(A4A5A6) 

90 perm 40 perm 

Y (AiA2A3)(A4A5)(A6)-i Y (AlA2A3)(A4)(A5)(A6) 

2 



— ^ z 

120 perm 40 perm 



+ Y (AlA2)(A3A4)(A5A6)-- E (AlA2)(A3A4)(A5)(A6) 

15 perm 45 perm 

+ i Y (AiA2)(A3)(A4)(A5)(A6) , (35) 

15 perm 

which is a useful relation for the construction of any 4-flavor chiral perturbation theory, such as [SU{A)l x SU{A)b. x 
U{l)v]/[SU{A)v X U{l)v], SU{A)/SO{A), and S't/(4)/S'p(4) chiral perturbation theories. 

The n = 5 case 

For any 5x5 matrix A, Cayley-Hamilton theorem reads 

= xa(A) = A5 - (A)A4 + i((A>2 - (a2))A-^ - i((A)3 - 3(A2)(A) + 2{A'))A^ 
1(A4) + i(A3)(A) + i(A2)2 - i(A^)(A)2 + ^(A)^) A 

\{A^) + i(A^)(A) + i(A-^)(A2) - \{A^){Af \{AY{A) 
^ 'A2)(A)3--i-(A)5') , (36) 



12' ' 120 

from which, the same procedure as the n = 2 case yields 

0= Y (A1A2A3A4A5A6) - Y (AiA2A3A4A5)(A6)- Y (AiA2A3A4)(A5A6) 

120 perm 144 perm 90 perm 

+ Y (AlA2A3A4)(A5)(A6) - Y (AlA2A3)(A4A5A6) 
90 perm 40 perm 

+ Y (AlA2A3)(A4A5)(A6)- Y (AlA2A3)(A4)(A5)(A6) 

120 perm 40 perm 

+ Y (AlA2)(A3A4)(A5A6) - Y (AlA2)(A3A4)(A5)(A6) 
15 perm 45 perm 

+ Y (^1^2) (A3) (A4) (As) (Ae) - (Ai) (A2) (A3) (A4) (As) (Ag) , (37) 

15 perm 

which is a useful relation for the construction of any 5-flavor chiral perturbation theory, such as [SU{5)l x SU{5)]i x 
U{l)v]/[SU{5)v X Uil)v] and SU{5)/SO{5) chiral perturbation theories. 
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IV. APPLICATIONS TO GENERAL RELATIVITY 



In general relativity, a second-rank covariant tensor field g^i, called the metric plays a crucial role. Its determinant, 
g = det((7pi,), is usually needed to be figured out in case one wants to handle the pseudoscalar volumn element: 
y^—gd'^x, where the superscript d stands for the number of spacetime dimensions. Let the dimension number d be 
10, which is required for the consistency of superstring theories. Using Eq. (jl5p . one can express the determinant 
of the metric 10 x 10 matrix g in terms of the traces of the successive powers of g. Since there are 42 solutions 
(pi,P2j • ■ • jPio) of the equation pi + 2p2 + • • • + lOpio = 10, there are also 42 terms in the expression of detg, which 
is given by 



det(g) 



1 

10 
1 



10 



-TgTi 

9 



1 



48 

'288 ^ 



24 



1 



i4- 
1 



l^TM - lr4T3 



1 



2880^ 

1 j,2j,6 1 



162 ^ 144 ^ 2 72 ^ ^ 



W^'^' ' W'^'^' ' 40^ 

1 

^ 432 



192 



5760 



80640 



15120 
1 

3628800 



96 



.10 



(38) 



where Tm {rn = 1,2,, 



, 10) stand for tr(g™). Eq. ([38|) shows that the determinant g = det(g) is completely decided 
by the 10 traces, Tm (m = 1,2, ... , 10), and this expression is fairly controllable and thus can be possibly used in 
future theoretical considerations. By using PASCAL (a sample code is given in App. [B|, one easily finds that there 
are 56 trace terms among det(g) for 11 dimensional spacetime corresponding to the M theory, while 2436 terms for 
26 dimensional spacetime. 



V. CONCLUSIONS 



In this paper, we have given a new proof, by using combinatorics, for the recurrence relations of the characteristic 
polynomial coefficients, and, moreover, we have obtained an explicit expression for the generic term of the coefficient 
sequence, which yields an improved version of the Cayley-Hamilton's theorem with all coefficients explicitly given, 
and which implies a byproduct expression for the determinant of any finite-dimensional matrix in terms of the traces 
of its successive powers. We have shown that these relations can be applied to chiral perturbation theories and general 
relativity. On the one hand, we have obtained the Cayley-Hamilton relations up to the trace of six operators in any 
n-flavor chiral perturbation theory for n runs from 2 to 5. On the other hand, we have obtained the determinant of 
the metric tensor in terms of the traces of its successive powers in 10 dimensional spacetime. For higher dimensional 
spacetime, it is also easily obtained by using a PASCAL program code, which is contained in the appendix. 



APPENDIX A: A CHECK FOR THE CAYLEY-HAMILTON RELATIONS IN THE n = 2 CASE 



In this appendix, we give a check for the Cayley-Hamilton relations in the n = 2 case by using the property of 
Pauli matrices. For dcfiniteness, let us verify Eq. (jSO]) in Example 2. Since is traceless due to Eq. ((27)) . then it 
can be decomposed, based on Pauli matrices, to = t'^X^ with the coefficients X^ (a ~ 1, 2, 3) and r" (a = 1, 2, 3) 
standing for Pauli matrices. Now let us show the following identity: 

{X^X.XpX,) = i (^{X^X,) {XpX,) + (X^X,) {X,Xp) - {X^Xp) {X,X„)^ , (Al) 

where (■ ■ • ) stands for taking the trace. 

Proof. The well-known identities t"t'' = 5"'' + ie°^'=r'= and e'''"' e"'^" = 5'"'S'"^ - S'"^6'"' imply that 
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Thus, 



{X^X^X.X,) = {T'^T\-r'')XlXlXlXt 



= 2{6''''S'"^ + S'"'6'"' - 5'"'5^'^)XlXlX'pX'^ 

= i (^{X^,X,) {XpX,) + {X^X,) {X,Xp) - {XpXp) {X,X„)^ . □ 
Then, by contractmg indices in Eq. (|Aip . we obtain 

{XpX^X^X'') = {X^X„Y - i(A^A^)2 , (A2) 
{XpX^X.X-) = ^{XpX'^f , (A3) 

from which it is easily checked that 

2{XpX'^X,X'') + {XpX.X^X-') = ^{XpX'^f + {XpX.f , (A4) 

which is exactly Eq. ([50)1 obtained in Example 2. It should be kept in mind that Eq. (|Aip . or Eqs. (jA2[) and 
(|A3p hold only for n = 2, though they arc stronger than Eq. (jA4p which comes from the Cay ley-Hamilton theorem. 
However, from Eq. ([5^ . one notes that the following relation holds as well in the n = 3 case. 

= (A1A2A3A4) - (AiA2)(A3A4) , (for (A,) = 0) . (A5) 

6 perm 3 perm 



As a consequence, Eq. (|A4[) holds not only for 2 x 2 but also for 3x3 traceless matrices X^. Obviously, this relation 
is invalid for n x n (n > 4) matrices. 



APPENDIX B: A PASCAL CODE FOR COMPUTATION OF D, 

Type splittype=array [1 . . 40] of integer; 
var num, sum, k: integer ; 
var spl : splittype ; 
var fl,f2:text; 
procedure Output; 

var i ,x ,y: integer ; 

var sp : splittype ; 

begin 

for i:=l to 40 do sp[i]:=0; 
for i:=l to k-2 do 
begin 

write (fl, spl [i] ,'+'); 
sp [spl [i] ] : =sp [spl [i] ] +1 ; 
end; 

writeCf l.spl [k-1] , ' = ' ,num) ;sp[spl [k-1]] : =sp [spl [k-1] ] +1 ; 
writeln(f 1) ; 

for i:=l to num do write (f 2 , sp [i] , ' '); 
i:=l; 

while sp[i]=0 do i:=i+l; 

if sp [i] 01 then write (f 2 , ' 1/ ' , sp [i] , ' ! ' , ' (-T_ i , V , i ,')"', sp [i] ) 

else write (f 2 , ' (-T_ ' , i , ' / ' , i , ' ) ' ) ; 
i:=i+l; 

while i<=num do 
begin 

if sp[i]<>0 then 

if sp[i]=l then write (f 2 , ' * (-T_ ' , i , ' / ' , i , ' ) ' ) 

else write (f 2 , ' *1/ ' , sp [i] , ' ! ' , ' (-T_ i ,'/', i ,')"', sp [i] * ') ; 
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i:=i+l; 
end; 
writelnCf 2) ; 
end; 

procedure split (n,m: integer) ; 
var i: integer; 
begin 

if n=0 then begin Output; suni:=sum+l end 
else if n>0 then 
begin 

for i:=m downto 1 do 
begin 

if n>=i then spl[k]:=i; 
k:=k+l; 
split (n-i , i) ; 
k:=k-l; 
end; 

end; 

end; 
begin 

sum:=0;k:=l; 

assignCf 1 , 'Output . txt ' ) ; 
assign(f2, '0utput2.txt') ; 
rewrite(f 1) ; 
rewrite (f 2) ; 
readln(num) ; 
split (num,num) ; 
writeln(fl , 'Total is ',sum); 
writelnCf 2, 'Total is ',sum); 
closeCf 1) ; 
close (f 2) ; 
end. 
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